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Forced transient motions of peripherally supported, circular, elastic 
plates are analyzed according to the classical plate theory. The Green's 
function for the plate is developed and used to construct solutions for con- 
centrated impulsive loadings, suddenly applied loadings, and moving 
pressure-wave loadings. The boundary of the plate is considered to be elas- 
tically built-in in a manner that prevents transverse edge motion and pro- 
vides a restoring edge moment linearly related to edge rotation. Thus, limiting 
cases include a clamped plate and a simply supported plate. Numerical 
results are included to illustrate the influence of structural and loading 
parameters on the dynamic response of the plate. 

I. INTRODUCTION 

During the past decade considerable effort has been channeled toward 
increasing the capability of equipment to sustain severe nuclear weapon 
environments. These efforts, commmonly called "hardening," employ 
various combinations of analytical and experimental approaches, each 
approach having certain difficulties and shortcomings. 

Problems that are analytically tractable are usually restricted to 
simple and regular geometries and usually incorporate simplifying ap- 
proximations as to material properties, weapon phenomenon, and sepa- 
ration of the combined weapon effects into independent and separate 
effects. One of the areas of hardening which lends itself to both analytical 
and experimental treatment and for which some full-scale nuclear test 
data are available is the response of structures to an air-blast wave; the 
analysis of structures subjected to air-blast pressure waves is the subject 
of the present article. In particular, the transient displacements of 
circular elastic plates subjected to impulsive and moving loads will be 
analyzed according to the classical plate theory. The classical, or small 
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deflection, plate theory does not account for internal structural damping, 
effects of transverse shear, or rotatory inertia. Consequently, this 
analysis is strictly applicable when the plate is "thin" (small thickness 
to radius ratio) and when the higher modes of vibration are of secondary 
importance;* neglecting internal structural damping results in predicted 
deflections that are conservative in the sense that they will be larger 
than corresponding deflections with damping present. 

Equations are derived and numerical results are presented for several 
elemental loadings and for pressure loading waves of constant and de- 
caying magnitude that sweep across the plate with uniform speed. Sweep- 
ing pressure-wave loadings occur when the blast wave approaches the 
plate in other than a face-on direction. 

The method of analysis is to construct the Green's function for the 
plate and then to use principles of superposition in synthesizing solutions 
for the various loadings of interest. 

Previous analyses of circular elastic plate vibrations can be traced 
to the free vibration analyses of Poisson, 1 and of Kirchhoff, 2 who con- 
sidered axisymmetric and nonaxisymmetric vibrations respectively. 
More recent studies of forced vibrations of circular elastic plates include 
the work of Flynn, 3 Sneddon, 4 Riesman, 6 and the present writer. 6 
Mindlin 7 discussed the effects of rotatory inertia and transverse shear 
deflections on the dynamic plate equation. Elastically restrained plates 
were investigated by Kantham 8 and by Reid. 9 

II. FORMULATION 

2.1 Equation of Motion 

Forced transverse motions of a homogeneous, isotropic, elastic plate 
of constant thickness (thickness is restricted to be small in comparison 
with radius) are governed by the partial differential equation 

D^wirM) + m d ^0- = V {r,6,t) (1) 

under the restriction that the deflections are small in comparison with 
the plate thickness, and that the shear deflections, rotatory inertia, and 
damping can be neglected. Fig. 1 defines the coordinate system and 
configuration that correspond to (1) 

* For a complete discussion of the limitations of the classical plate equation, 
the reader is referred to Ref . 7. 
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Fig. 1 — Circular plate configuration. 



2.2 Boundary and Initial Conditions 

A plate that is continuous at the origin (non-annular) and elastically 
built-in* along its periphery has the four boundary (or regularity) 
conditions: 



(a) 
(6) 
(c) 

(d) 



w(0,d,t) must be finite 

dw(0,d,t) , B . 

must be finite 

dr 

w(a,6,t) = 

d~w(a,d,t) 



-D 



dr- 



v dw{ 

a 



[a,6,t) ~\ 

dr J 



= fi 



dw(a,6,t) 
dr 



(2) 



Initial conditions that characterize an initially undeflected and sta- 
tionary plate are 



w(r,0,O) = 



and 



dw(rfifi) 
di 



= 0. 



(3) 



These initial conditions are sufficiently general for the applications to 



* Elastically built-in edge, as used in this article, refers to a boundary support 
that prevents transverse edge motion and provides a restoring edge-moment pro- 
portional to edge rotation. By properly selecting the constant of proportionality, 
special cases corresponding to a clamped edge and a simply-supported edge are ob- 
tained. A more detailed discussion of this boundary condition appears in Ref. 6. 
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be considered here; other initial conditions can be treated in a similar 
manner and will introduce minor changes in the equations developed 
herein. 

III. BASIC SOLUTION 

3.1 Dimensionless System of Equations 

Introducing the dimensionless quantities p = r/a, i) = w/a, r = t/T, 
a = ma*/DT 2 , and 7 = v + (0a/ D), (1) and (2)d become, respec- 
tively, 

V i)(pfi,T) -(- a — - = -=r p{p,9,T) (4) 

a-r D 

and 

aVur) ., jjOjgj _ a (5) 

dp 2 dp 

3.2 Homogeneous Equation 

As a prelude to the solution of (4), the set of eigenf unctions for the 
plate will be determined. These eigenfunctions arc obtained by starting 
with the homogeneous equation 

vVpAr) + a *&£> = 0. (C) 

or* 

Separable product solutions of such form that the angular functions 
and time functions possess the necessary periodicity suggest the follow- 
ing functions: 

!cos n9\ 
or K U "' T , (7) 

sin nd\ 

where i = s/— 1 , and n and j are integers. 

Substituting (7) into (6) produces an ordinary differential equation 
in p 

(£+U-?+^)(£ + ;i-? -•')*-«-* (8) 



where 



2 i 

K n j = a 0)„j 
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Equation (8) and the boundary conditions (2) define radial eigenf unc- 
tions for the plate. 

These radial eigenfunctions are readily found to be 

R nj ( P ) = In(Knj)Jn(KniP) " /.(**)/.(«■#), (9) 

where k„j is determined from the transcendental equation 

^- InUni)Jn(lCni) = h U nj )J n+l U nj ) + S,(«rf)W*i)- (10) 
1 - 7 

Two identifying subscripts are associated with k„j because there is a 
doubly-infinite set of eigenvalues and eigenfunctions. The first subscript, 
n, indicates the number of nodal diameters occurring in that mode of 
vibration and can be any positive integer (representing nonaxisy in- 
metric modes) or zero (representing axisymmetric modes); the second 
subscript, j, indicates the number of nodal circles (including the bound- 
ary circle) and can be any positive integer. Eigenvalues for representa- 
tive parameters are presented in Table I. 

The radial eigenfunction R n j(p) as defined by (9) and (10) will be 
used in the next section as a building block for the development of the 
Green's function. 

3.3 Green's Function 

Next, the Green's function associated with (4) will be developed. 
Toward this end, consider the loading function 

p(p,0,r) =- 8(p - p )8(d)8(r - r„). (11) 

Po 

The solution of (4) with the particular loading function (11) is de- 
fined as the Green's function and is denoted by G(p,6,t; p a , 0, t ). With- 
out loss of generality, G may be represented by the double series 

00 00 

Q(p,6,t; po , 0, to) = E E Rv(p) cos kd g kj {r), (12) 

where 0*>(t) is an unknown function of time that must now be deter- 
mined. 

Introducing (12) into (4), and interchanging the order of differentia- 
tion and summation, the result is 

E E if, + - T ~ tY 8*00 cos ke g kj (r) 
a- , W p dp p 2 / 

Y" V / » I \ ia ■■ i \ a hip — Po)S(6)8(t — To) 

+ «LL Rhi(p) cos kO g kJ (T) = jz . 

k j U Po 
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Capitalizing on the fact that Rkj(p) satisfies (8), the first differential 
operator in (13) can be replaced by **/ and (13) becomes 

£ £ k-;W r ) + <xg kj (T)]R ki ( P ) cos kd 

k 

_ a 3 5(p - p )8(d)8(r - r,) (14) 
D po 

Next, both sides of (14) are multiplied by 

pRim(p) cos Id dddp, 

and integration is performed, first with respect to from — > 2ir and 
then with respect to p from — > 1 ; by virtue of the orthogonalities (see 
Appendix), this reduces to 

3 

nfruNrnJlUr) + aQnNvJgiM = ^R lm ( P o)8(r - r ), (15) 
where 

7T (I 9* 0), 



011 - \2tt (Z = 0), 
and 

Nmm = h[h (>Clm)Jl+l (*lm) ~ Jl («Jm)/j+l («Im) } 

fl +7 + 2Z1 -2/ w2, v 

— : ll {Kl m )Jl (Kim). 

. 1-7 J 
For the prescribed initial conditions, the solution of (15) is 

, v a 1 Rkj(po) SU1Mw(t — To) , s (AR s 

QkM = n - r, v k Kt - To). U6J 

Equation (16), when used in conjunction with (12) is the Green's 
function for the circular elastic plate. Physically, it is the response due 
to a transverse point impulse applied to the plate. 

3.4 Generalized Green's Function 

The Green's function given by (12) was developed for the case of an 
impulsive loading singularity located on the line 6 = 0. Later, need will 
arise for the response due to a loading singularity located at the point 
(p , — a Q ) (see Fig. 2). Employing the trigonometric identity, 

cos (.4 4- B) = cos A cos B — sin A sin B, 
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Fig. 2 — Plate coordinate system. 

a trivial change of variables produces the corresponding generalized 
Green's function G(p,8,t; p a , a , t ) 

00 00 

G(p,6,t; po , «o , To) = EE Rkj(p) cos ka a cos he </*>(t) (17) 

00 00 

- X) 2 Rkj(p) sin ka sin /e0 </*,(r). 

A=l |-1 
IV. APPLICATION OF GREEN'S FUNCTION 

The Green's function as given by (17) will now be utilized to construct 
solutions for several loadings of technical interest and of practical con- 
cern. 

4.1 Ring Loading 

The first example is that of an impulsive concentrated ring pressure 
distribution given by 

v(pAt) = pW 6( " - "° )S{T ~ To) (18) 

Po 

(see Fig. 3). Physically, this loading corresponds to a concentrated ring 
impulsive loading applied at radius p and at time t . A differential 
element on the ring has a force per unit length p(p )/p . The resultant 
displacement of the plate is obtained by superposing contributions due 
to all of the elements on the ring. This superposition is expressed mathe- 
matically by the integral 

v(p,0,t) = I G(p,6,t; p ,(Xo,To)p(po)da 0i (19) 

Jo 

where a is defined by Fig. 3. 



TRANSIENT MOTION OF CIRCULAR ELASTIC PLATES 2417 




Fig. 3 — Ring loading. 

The evaluation of (19) may be simplified considerably by observing 
that for I 9^ the integrals are all zero, leaving only the I = terms. 
Accordingly, (19) simplifies to 

( v ap(po) v^ Roj(p)Roj(po) . / \,/ \ / on \ 

v(p,d,r) = — ?-=- 2-, TT~ « sinaioi(T — t )1(t — t ). {20) 

Da pi Njfuoj 

Equation (20) may be recognized as the axisymmetric solution that 
was obtained previously, 6 provided of course that proper account is 
taken for the change in nomenclature. 

4.2 Concentrated Line Loading 

As a second example, consider a concentrated line impulse loading 
specified by 

P (x,y,r) = S(x- X )8(r - r ), (21) 

where X and y are cartesian coordinates appropriate for this problem and 
defined by Fig. 4. Physically, this loading corresponds to an impulsive 
concentrated line loading applied at position x = x a and at time t = t„ . 
The force per unit length acting on each element of the line loading 
given by (21) is 

«(r - T ). 

Likewise, the force acting on the plate due to the Green's function 
singularity located at (.t„ , y) is 

5(T - T ). 

It follows that the resultant deflection due to all differential elements 
of the line loading is expressed by the superposition integral 
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y 




Fig. 4 — Concentrated line loading. 



/.(l-x 2)» / \ 

r,(p,d,r; Xo , T ) = / G( P ,d,r; [x 2 + y]\ tan ' M- , T ) dy. (22) 

Substituting (17) into (22) gives the following expression for the 
dynamic response of the plate: 



j?(p,0,t; x ,To) 



•(l-*o 2 )* 



Da J~0 m=l L WllNmm'oOlm •> -(l-i 2)i 

r, , -i y~\ , fl*m(p) sin Z0 [ {l ~ x ° ) D , 2 , „ «v* 
•cos Z tan —dy ^— £ — ; / Rim{x +y ) 

[_ X J WiiN mm l O)i m J-(i- Io 2)J 

• sin I tan -1 - d^ X sii 



(23) 



Sin U3 lm 



(t — r )l(r — To). 



The integrals in (23) are not expressible in closed form; however, 
they can be integrated numerically. Fortunately, there is one point on 
the plate where the evaluation does simplify considerably, and that is at 
the center of the plate, i.e., p = 0. Restricting our attention to the center 
of the plate, only the I = terms are nonzero, and consequently the 
center deflection is given by 



„(0,t;.t o , To ) =^Z n 



I 

Jet 



O30m J 

X sin uomir — r )l(r — r Q ). 



Romixo 2 + 2/ 2 ) 1 dy 



(24) 



The integrals in (24) are perfectly well behaved and finite, although 
they cannot be expressed in closed form. Certain special cases (e.g., for 
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20 30 40 

TIME IN MILLISECONDS 



Fig. 5 — Dynamic response due to impulsive line loading (y — 1.0). 

Xo = 0) have been tabulated. 4 Several other cases have been evaluated 
by the writer, and the results have been incorporated into the numerical 
examples. 

4.3 Pressure Acting on a Portion of the Plate 

A third example is a plate that is impulsively loaded by a uniform 
pressure acting on a segment of its surface and unloaded elsewhere. 
This loading is depicted by Fig. 7. 

The center deflection may be thought of as the resultant deflection 
due to loadings on each line segment from 

— 1 ^ Xo ^ x . 

A pressure of magnitude P Q S(t — t ) loads each element of width dx a 
with a line loading of magnitude P 8(t — t q ) dx . The resultant center 
deflection is obtained directly from the superposition integral 
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Fig. 6 — Dynamic response due to impulsive line loading (7 = 0.33). 




Fig. 7 — Partially loaded plate. 
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t)(0,t; x',t ) = / P o y(0,T; Xo,t ) dxo, (25) 

where tj(0,t; x a , r D ) is defined in (24). 

The integrals in (25) may be integrated numerically. One case that 
is readily integrable in closed form is when x = 1 ; that case corresponds 
to a plate that is suddenly loaded by a uniform pressure over its entire 
surface. It is readily shown in that instance that (25) reduces to the 
previously obtained solutions presented in Ref. 6. 

If a sector of the plate is loaded with a uniform pressure impulse load- 
ing, as indicated by Fig. 8, the superposition methods employed above 
confirm the intuitive suspicion that the center deflection is the deflection 
for a uniformly loaded plate multiplied by the quantity 6 /2ir, where 
d a is the included angle of the sector (expressed in radians). 

4.4 Sweeping Pressure Wave 

Next, consider the case of a circular plate loaded by a step blast wave 
that sweeps across its surface with constant velocity c (expressed in 
plate diameters per dedimensionalized time) as indicated in Fig. 9. The 
sector of the plate behind the wave front is loaded by pressure P Q , 
whereas that portion ahead of the wave front is as yet unloaded. 

It is convenient to consider the plate to be divided into strips of equal 
width A.r , with the strips numbered successively from the left hand 
edge (x = —1)- This is depicted in Fig. 10. 

The first strip is then loaded by a step pressure wave (unloaded until 
T = 0) and then loaded with pressure P for r D > 0). After an increment 




Fig. 8 — Plate loaded on sector. 



2422 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1965 




Fig. 9 — Sweeping pressure wave with velocity c. 

of time At has elapsed, the second strip is loaded by the same step 
pressure wave (unloaded until t d = At„ and then loaded with pressure 
Po for To > At ). After time interval 2At c , the next strip is loaded by 
the same step pressure loading, and so on, until eventually all strips are 
loaded by pressure P . As the width of the strip is reduced to an in- 
finitesimal, and correspondingly the number of strips increases to in- 
finity, the continuously sweeping pressure wave results. 

The solution for a suddenly applied step line loading is obtained by 
integrating (23) with respect to t g (Duhamel integral). The suddenly 
applied step line loading, denoted by v(p,6,t; x , t ), is given by 

7j(p,d,T;x ,To) = / r}(p,d,T; x , To) dr . 




Fig. 10 — Sweeping pressure wave. 
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The result of this integration is 
rj(p,d,T;x ,To) 



a ^ ^ Rimip) cos 16 f p / i , • 

xA dy ~ e H AT-W i- 



cos Z tan — dy — 



S2m(.-Co + y ) 



§ (26) 



sin Ztan -1 -^- LJx[l- 



ttjro* J-(l-x *)l 
COS co i m (r — t )]1(t — To). 



In terms of the solution for the suddenly applied step line loading, the 
sweeping step wave solution, as obtained by superposition is 



V*(p,6,t; c) = lim J-, PoV ( pfi,r\ - 1 + nAx , - 

Ax n -*0 n=0 \ C 



! W , 



(27) 



for 0^t|- 
c 



where ij is given by (26), and rj* is the deflection due to the sweeping 
wave. In the limit, (27) is expressed by the integral 



r,*(p,d,r; c) = P f v (p,d,r; y , L±JM dy ^ r ^ -. 
Jy u =-l \ C / C 



(28) 



Although (28) cannot be evaluated in closed form, (27) is perfectly 
well suited for numerical evaluation (see Figs. 11, 12, 13), and the ac- 
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Fig. 11 — Dynamic response due to sweeping step wave (7 = 1.0). 
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0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 2.25 2.50 

PLATE DIAMETERS TRAVELED BY BLAST WAVE FRONT 

Fig. 12 — Dynamic response due to sweeping step wave (7 = 0.33). 

curacy obtainable is limited only by the computer time expenses associ- 
ated with the numerical integrations. In the numerical evaluations 
included in this article, certain symmetry properties were exploited to 
minimize computer time requirements. 

4.5 Sweeping Pressure Wave with Decay 

As a final example, the case of a sweeping pressure pulse will be con- 
sidered where the magnitude of the pressure at any point behind the 
wavefront decays exponentially as a function of the distance behind the 
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Fig. 13 — Dynamic response due to sweeping step wave (y = ■»). 
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wavefront. As in the previous example, the speed of the pressure wave is 
considered to be constant. 

Referring to Fig. 14, a point at x = x Q is unloaded until time 
T = (l -f Xo )/c, at which time the pulse arrives, loading it with pres- 
sure P . Because of the decay of the loading pulse, the point at x is 
loaded with pressure 



P(So,r) - P exp {| in |[ T - L±^]} 1 (r - L±^) . 



(29) 



The expressions for the deflection due to the loading given by (29) are 
identical to (27) and (28) with the exception that rj must be modified. 
The necessary modification is to replace the factor [1 — cos oj im (T — r G )] 
that appears in (20) by the corresponding factor 



O-'/r 



In 5 



,<!exp [1ii5(t -- t )] - sin wi«(t — T ) 

coi„, 2 4- (In 6)- { wim 

— COS U31„,(t — T )>. 



(30) 



It can be seen that for the case of zero decay (In 5 = 0) (30) reduces 
to the step loading solution (Figs. 15, 16, 17). 

V. RESULTS 

Dynamic response curves are presented in this article for line impulse 
loadings and for sweeping waves of both constant and time-diminishing 
pressure pulses. Figs. 5 and 6 illustrate that a line impulsive loading ap- 
plied near the diameter excites primarily the first mode of vibration; line 
loadings applied away from the diameter excite a larger percentage of 
the higher modes, although as one should expect, the magnitude of the 
deflection diminishes as the line loading is applied nearer to the edge. 




Fig. 14 — Sweeping pressure wave with decay. 
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functions were computed. In this appendix, corresponding orthogonality 
properties and normalizing factors will be explored for the non-axisym- 
metric modes of vibration. 
The eigenfunctions, 

[cos nd) [cos nd) 

R„j(p)\ Or \ = {I n (Knj)Jn(Knjp) ~ Jn(Knj)In(K nj p)} \ OT k 

[sin nd) [sin nd) 

where n = 0,1,2, • • • and K nj satisfies the transcendental equation 






In(K„j)J„(K n j) = I n (K„j)J n ^l(K„j) + Jn(Knj)In+l(l<nj), 



are solutions to the partial differential equation (6) and satisfy the 
boundary conditions. 

A.l Orthogonality 

To begin, recall the well known 10 results 



f 



cos nd sin m9dd = 



cos nd cos mddd = < 



m,n any integers; 



r 


(n ^ m) y 


7T 


(n = m t* 0), 


[2lT 


(n = m = 0); 


[° 


(n 7* m), 


. 


(n = m). 



I sin nd sin mddd = 
Jo 

These orthogonality relationships imply, in our case, that there is no 
coupling between modes of vibration having n nodal diameters and those 
modes having m nodal diameters, where n t± m. Thus, one need only 
concern oneself with relationships such as 



/ pR n] (K„jp)Rnk(ii n kp) dp. 

Jn 



Therefore, consider the ordinary differential equations (8) that R nj - 
and Rnk must satisfy 



/ d , 1 d n\ p 4 P 

\dp 2 pdp p 2 / 



and 
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{h + \i~v) R " k = KjRi "" 



Multiplying the first equation by pR nk dp, and the second by pR„jdp, 
subtracting, and integrating with respect to p from — > 1 gives 

(«!./ — Knk) I pRnjRnkdp = I pR„k [ j~S + ~ j o ) ^»3^P 

Jo Jo \dp 2 prfp p 2 / 

The right side of (32) is expanded, and integrated by parts. The 
lengthy, but routine details of this operation are suppressed in the inter- 
est of brevity; the result is 



/ pR, t jlt„kdp 

Jo 



f lPR Hi (PRnj 1 dR„i\ j tPRnk (PRnk 1 dRnh 

Bnt \ p 1 P r + ~d^~'p dp~)~ Bni Y d P > + dp* ~ P d P 

[ dRnid*R«h dRnk&RnA ,2n a jdRnk n dR ni „ 

X l P dp dp* dp dp 2 J p \ dp dp 



(33) 



W - K„k*) 

For the case of n = 0,(33) was shown to vanish in Ref. G forj ^ k. For 
n = 1,2, • • • similar reasoning demonstrates that each of the terms in 
(33 ) either vanishes or the terms mutually annihilate each other, so that 
the eigenfunctions are orthogonal under the rather general boundary 
conditions given in (2). 

A.2 Normalizing Factor 

Forn = 0,1,2, • • • , but where j = k, both numerator and denominator 
of (33) are zero. Thus, a limiting process must be employed to find the 
value of 



Njj" m ( pR„/(p)dp 



from (33). Differentiating numerator and denominator of (33) with 
respect to k„* and then setting K nk = k»j , the result is 

Njj" = h{In(Knj)Jn+l(Knj) — Jn U„j)I n+ l (k„j)} 

1 + t + 2n T „ , T 2 ( s (34) 

: in UnjVn \K n j) • 

1 - 7 
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This normalizing factor is identical to the axisymmetric normalizing 
factor derived in Ref. 6 for the axisymmetric case (n = 0). An alterna- 
tive way to derive (34), which was used to check this equation, is to 
perform the indicated integration directly. 

NOMENCLATURE 

a radius of boundaiy of plate 

c speed at which pressure wave sweeps across plate (ex- 

pressed in plate diameters per dimensionless time unit) 

(Eh 3 \ 
= 19 /i _ o\ ) 

E Young's modulus 

G Green's function 

gkj(r) function of time defined by (10) 

h plate thickness 

Jn(z),I n (z) nth order Bessel function and modified Bessel function of 

argument z 
m mass per unit projected area of plate 

JV/y" normalizing factor for radial functions, defined by (34) 

p(pA T ) pressure loading on plate 

P pressure of loading pulse (constant) 

Rnj(p) radial eigenf unction defined by (7) 

r radial coordinate 

t time 

T time duration used to de-dimensionalize the time variables 

(chosen for convenience) 
w plate deflection 

x,y cartesian coordinates 

dimensionless parameter (= ma /DT ) 



a 



/3 modulus of spring on edge of plate 



edge-fixity parameter ( = v + — J 



Si 

\ ' ' D 



8 decay rate of sweeping pressure pulse 

8(z) Dirac delta function of argument z 

■q dimensionless plate deflection (= w/a) 

8 angular coordinate 

®n normalizing factor for angular normal functions 

K n j eigenvalue corresponding to mode with n nodal diameters 
and j nodal circles 
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v Poisson's ratio 

p dimensionless radial coordinate (= r/a) 

t dimensionless time variable (= t/T) 

u n j dimensionless angular frequency associated with the n—j 

mode of vibration 

V 2 Laplacian operator (= ±- + 1 ± + ± ^) 

1 (z) unit step function of argument z 

Differentiation with respect to time is denoted by dots. 
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